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Abstract: We discuss our model for soft colour interactions and present results 
on the diffractive structure function F®(f3,xjp,Q 2 ) and inclusive transverse energy 
. flows, which agree with available HERA data. 
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The observation of rapidity gap events at HERA Jl|, || has led to a renewed interest in 
diffractive phenomena and how they can be understood within QCD. In || we presented a 
new model for diffractive hard scattering based on a mechanism for soft colour interactions 
(SCI). The starting point is the normal DIS parton level interactions, with perturbative QCD 
corrections based on standard first order matrix elements and conventional leading log parton 
showers for higher order parton emission. This gives a state of colour-ordered partons that 
will have further non-perturbative interactions to produce the final hadron state. The novel 
feature in our model is the introduction of random soft colour interactions corresponding to 
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non-perturbative gluon exchange between these partons. This changes the colour structure and 
thereby affects the hadronic final state when a conventional hadronisation model such as the 
Lund string model 0] is applied. 



> 



Qh! Rapidity gaps may then arise when a gap at the parton level is not spanned by a string. In 
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particular, when the hard process starts with a gluon from the proton, leaving a colour octet 
remnant and a colour octet hard scattering system, a soft colour exchange between the two octet 
systems can give two colour singlets separated in rapidity ||. Large forward rapidity gaps are 
here favoured by the large momentum of the remnant, in particular at small-x. Our model || 
is implemented in the Monte Carlo (MC) Lepto 6.4 producing complete events. 

To compare our model with experimental data on rapidity gap events we consider the diffrac- 
tive structure function F-f xjp, Q 2 ) defined by [|| 
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(i.e. assuming single photon exchange and neglecting Fl). This inclusive quantity contains 
the dependence on the main variables f3 ~ Q 2 /(Q 2 + M x ), xp ~ (Q 2 + M X )/(Q 2 + W 2 ) and 
the usual DIS momentum transfer squared Q 2 . The acceptance corrected data from HI P] are 
compared in Fig. 1 to our model results obtained by selecting MC events with rapidity gaps 
similar to the HI definition (i.e. no energy in r\ max < rj < 6.6 where rj max < 3.2). 

The model is generally in good agreement with the data. It has a tendency to be below the 
data at large Q 2 , possibly due to slightly too much parton radiation. The /3-dependence seems 
to be the same in the model and the data and thereby the Mx dependence is also basically 
correct. The ayp-dependence in the model may be slightly steeper than in the data. Fitting a 
universal :r^, a -dependence we get a = 1.5 ± 0.2 from our model to be compared with the HI 

1 Presented by G. Ingelman at Workshop on 'HERA Physics - Proton, Photon and Pomeron Structure', 
Durham, September 1995, to appear in the proceedings. 
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result a = 1.2 ± 0.1 § and ZEUS results a = 1.3 ± 0.1 |] and a = 1.5 ± 0.1 § using different 
methods. However the result obtained in the model depends significantly on which xjp-bins that 
are included in the fitting procedure. This originates from the model not quite giving a straight 
line in this plot, but having a tendency for a curvature with smaller slope at small xjp and larger 
slope at large xjp. With improved experimental data this aspect of the model can be tested. 
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Figure 1: The diffractive structure function -F 2 D (/3, xjp, Q 2 ) from our soft colour interaction 
model (histograms) compared to the HI data points 0]. 



Our model is meant to be applicable for DIS in general and should therefore also be compared 
with normal DIS events. The observed large forward transverse energy flow in an inclusive 
event sample requires a substantial energy and particle production and is thereby 'orthogonal' 
to forward rapidity gaps. In Fig. ^ we compare our model result with the HI data ||. The 
agreement is quite good except for the smallest x-values where the model is below the data in 
the central region of the hadronic cms. 
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The possibility to obtain rapidity gaps depends on to which extent there are gaps at the 
parton level, i.e. on the amount of parton emission in the forward region. In our model we 
use the first order QCD matrix elements for the primary emission and therefore the treatment 
of their soft and collinear divergences is of importance. We have replaced the conventional 
requirement mfj > y cu tW 2 on any pair ij of partons with an advantageous 'mixed' scheme using 
cuts in s = (pi + P2) 2 and z = p\P/qP which handles initial state collinear singularities better 
||. Numerically we use s m - m = 1 GeV 2 and z m j n = 0.01. 
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Figure 2: Transverse energy flow versus pseudorapidity rf in the hadronic cms from our model 
(histograms) compared to HI data points || . 



Higher order parton emissions are taken into account approximately through parton showers 
developing from the incoming and scattered parton. The initial state shower is most important 
for the forward rapidity region. In our model we are using the conventional GLAP evolution 
scheme [10] summing leading logQ 2 terms. Terms with log(l/a;) are neglected in GLAP, but 
will become important at small x and should then be resummed as in the BFKL equation Q. 
Therefore GLAP evolution should no longer be valid at some small x. Recent studies [^] based 
on the CCFM equation which sums both leading logQ 2 and log(l/x) terms suggests that 
this happens only at very small x, below the region 10~ 4 ^ x ^10 -2 where the rapidity gap events 
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are observed. Thus, there is no strong reason that the GLAP formalism should not be applicable 
for our purposes. 

One may worry that unsuppressed parton emission will destroy the gaps ||. The cut-off in 
parton virtuality, defining the borderline to the non-perturbative region, is a regulator of this. 
We have previously @ used the value 2 GeV, but in our improved model || returned to 1 GeV, 
which has been standard and conforms with e.g. LEP analyses. Thus, our model works with a 
conventional unsuppressed GLAP parton shower. 

Of importance for the gap rate is also the fluctuations in the initial parton emission |J. 
Although, one may expect that a GLAP parton shower gives a fair mean description of events, 
there is no guarantee that it accounts properly for fluctuations. Larger fluctuations of the 
number of emitted gluons would increase the rate of gap events and also increase the inclusive 
forward energy flow due to 'downwards' and 'upwards' fluctuations, respectively. 

In summary, our model based on non-perturbative soft colour interactions can explain im- 
portant features of both rapidity gap events and normal DIS interactions at HERA. 
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